On thermalization in Kitaev's 2D model 



R. Alicki^'^ M. Fannes* and M. Horodecki^'" 

^ National Quantum Information Centre of Gdansk, Poland 

^ Institute of Theoretical Physics and Astrophysics 
University of Gdansk, Poland 

^ Instituut voor Theoretische Fysica 
K.U.Leuven, Belgium 

Abstract 

The thermalization process of the 2D Kitaev model is studied within the 
Markovian weak coupling approximation. It is shown that its largest relax- 
ation time is bounded from above by a constant independent of the system 
size and proportional to exp(2A/A;T) where A is an energy gap over the 4- 
fold degenerate ground state. This means that the 2D Kitaev model is not 
an example of a memory, neither quantum nor classical. 

1 Introduction 

The fragility of genuine quantum states in the presence of interaction with 
an environment, in contrast to those which allow a classical interpretation, 
represents the main challenge for the large scale implementation of the ideas 
of quantum information processing. 

It is well-known that to improve the stability of quantum information 
processing the logical qubits should be implemented in many-particle sys- 
tems, typically N physical spins per logical qubit. The logical qubits should 
be stable objects with efficient methods of state preparation, measurements 
and application of gates. By efficiency we mean certain scaling behaviour, 
e.g. the lifetime of a logical qubit should grow exponentially with while 
the number of steps needed to apply a gate and a measurement should at 
most increase polynomially with A^. 

There exist different sources of noise, some of them of macroscopic ori- 
gin and often highly system-dependent. They can, at least in principle, be 



substantially reduced by proper engineering. However, one source of noise is 
inescapable: the microscopic interactions of the physical spins with thermal 
particles or excitations of the local environment. 

The analogous situation for classical information processing is well-under- 
stood. E.g. ferromagnetic interactions stabilize bits encoded is small mag- 
netic domains which can then be flipped by external magnetic fields. The 
existence of a similar mechanism for quantum information is still open to 
debate. 

As a first step, one is interested in designing a stable quantum mem- 
ory, i.e. a A^- particle system which can support at least a single encoded 
logical qubit for a long time, preferably growing exponentially with A^. The 
most promising theoretical candidates are provided by quantum spin systems 
which exhibit topological order, such as the 2D [1] and 4D Kitaev models [2] . 
The issue of thermal stability of these models has recently attracted much 
attention (see e.g. [3, 4, 5, 6]). The aim of this paper is to rigorously analyze 
the thermal stability of the 2D Kitaev model. 

The quantum theory of open systems developed more than thirty years 
ago provides a natural framework for studying stability in the presence of 
thermal noise. The particularly simple properties of Kitaev's model allow 
to apply the most developed tools in the theory of quantum open systems 
namely Davies's theory of semi-groups of completely positive maps which 
describes the dynamics of a quantum system weakly interacting with a heat 
bath in the Markovian approximation. Such evolutions satisfy the quantum 
detailed balance condition which permits a rather detailed analysis of the 
relaxation properties of the relevant degrees of freedom. 

We first develop some technical tools which provide upper bounds for 
the relaxation times of observables. Then we apply this to the study two 
examples: a ID quantum Ising model, which can be seen as a ID Kitaev 
model, and the 2D Kitaev model, both coupled to a heat bath. We show 
that the relaxation times of the observables which are supposed to encode 
the qubits are bounded from above by a constant, independent of the system 
size. Therefore, neither the ID nor the 2D model provides a stable quantum 
memory. Although this result was expected, the most important issue is the 
mathematical technique which is developed here, it could be applied to ana- 
lyze other models such as the 4D Kitaev model, as well. The method relies 
on a Hamiltonian picture of the dissipation where estimating the longest re- 
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laxation time is shown to be equivalent to estimate the ground state spectral 
gap for a very particular quantum spin Hamiltonian. These techniques, to- 
gether with those developed in our previous paper [4] on this subject could 
be applied to analyze other models. 

The paper is organized as follows: in Section 2 we briefly remind the ori- 
gins and the main properties of Davies generators as the quantum analogues 
of classical Glauber dynamics. Section 3 relates the lifetime of information 
encoded in a relaxing system to the spectral gap of the dissipative generator. 
In Section 4 a Hamiltonian picture of relaxation is obtained. The generator 
is recast into a master Hamiltonian on Liouville space such that the gap of 
the generator is recovered as the ground state gap of the master Hamiltonian. 
Section 5 deals with a simple model: the ID Ising model on a ring which 
can encode 1 topological qubit in its ground state. A rigorous argument is 
developed to show that the stability of this model when coupled to a thermal 
environment does not increase with the system size. Finally, in Section 6 
Kitaev's 2D model is shown to behave in a similar way as the ID Ising model 
with respect to thermal stability. 

2 Davies generators 

Davies generators provide a simple and realistic description of a system 
weakly coupled to a thermal environment. The interaction with the out- 
side world is so small that the system doesn't loose its identity and can still 
be considered on itself. The footprint of the thermal environment is in the 
reduced dynamics: a Markovian semi-group of unity preserving completely 
positive maps of a very specific type with the temperature encoded. The 
aim of this section is to briefly recall the general form and properties of such 
maps, henceforth called Davies maps, see [7] for the original rigorous analysis 
of the weak coupling limit. 

A small system, with a flnite number of states, is coupled to one or more 
heat baths at a same inverse temperature (5 leading to a total Hamiltonian 

H = H^^^ + i/^^*'^ + with = ^S^® fa, (1) 

a 

where the Sa are system operators and the bath operators. Both the 
coupling operators and are assumed to be Hermitian. An important 
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ingredient is the Fourier transform g^i of the auto-correlation function of 
/„,. The function ^„ describes the rate at which the couphng is able to 
transfer energy between the bath and the system. Often a minimal coupling 
to the bath is chosen, minimal in the sense that the interaction part of the 
Hamiltonian is as simple as possible but still addresses all energy levels of 
the system Hamiltonian in order to yield an ergodic reduced dynamics. The 
necessary and sufficient conditions for ergodicity have been obtained in [9, 10] 

{S^,H'^']' = £1, (2) 

i.e. no system operator apart from the multiples of the identity commutes 
with all the ^„ and H^'^''. 

The weak coupling limit results in a Markovian evolution for the system 
'^-^ = g{X)=i5{X)+L{X). (3) 

The generator is a sum of two terms, the first is a usual Liouville-von Neu- 
mann derivation as in standard quantum mechanics while the second is a 
particular type of Lindblad generator 

5{X)^[H'^\X] (4) 

. . ^ r . ^^t r.. ^ / M . r/^ / ^^t ,.1 ^ . X 

iuj) 



e-^- S^{u) [X , {S^{u))^] + e-^- [S^{u) , X] {S^{u))'^]. (6) 



+ 



Here the Saioj) are the Fourier components of as it evolves under the 
Hamiltonian system evolution 



^ae-*^^" = J]5„(a;)e--*, (7) 



where the a;'s are the Bohr frequencies of the system Hamiltonian. The 
Hermiticity of is equivalent with 

{S^{uj))^ = S^{-uj). (8) 
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Clearly, the temperature of the environment appears in (6) through the Boltz- 
mann factor. Such generators are called Davies generators in the sequel. 

A super-operator £ as in (3) generates a semi-group of completely posi- 
tive identity preserving transformations of the system. However, due to its 
specific form, it enjoys a number of important additional properties 

• the canonical Gibbs state is stationary 

Tr U e*^(X)) = T,{pp X) with pp = — ^, (9) 

• the semi-group is relaxing: any initial state p evolves to pp 

lim Trfpe*^(X)') = Tr(pflX), (10) 

• each of the terms Cauj, and therefore JC as well, satisfies the detailed balance 
condition, often called reversibility 

[S,£a^] = and Tr(^pf3Y^ LUX)) = Tr(p^ {C^Y))^ X^ (11) 

Equation (11) expresses the self-adjointness of Cauj with respect to the 
Liouville scalar product 

{X,Y)fs:=TrppX^Y. (12) 

Writing Q = iS + C, sec (3), is therefore decomposing C into a Hermitian and 
a skew-Hcrmitian part and, by (11), Q is normal. The dissipative terms C^uj 
of the generator, being self-adjoint, are negative definite. This implies that 
we can always sandwich the spectrum of C between spectra of generators 
where the g{u!) have been replaced by constants. 

The negativity of Cau; becomes manifest if we write 



(13) 



-{X , CUX))^ ={ [Saico) , X] , [S^iu) , X] >^ 

+ e-^-(K(^),X], [Si{u;),X])^. 

To show this equahty one uses that 

pp Saiio) = c'^ Sa{uj) pp, (14) 
which is a direct consequence of (7). 
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3 Relaxation times and memory 



Due to ergodicity any initial state of a system whose dynamics is governed 
by a Davies generator will eventually relax to equilibrium. Information can 
be encoded by perturbing the equilibrium state of the system and, in order 
to retrieve this information, one has to single out observables that detect the 
perturbation of the state. The longer lived these observables, the more effi- 
cient the memory. The aim is indeed to increase the useful life of observables 
and states by devising proper Hamiltonians and encoding procedures. We 
are in this paper interested in the 2D Kitaev model of size N and especially 
in the relation between lifetime and system size. We begin by arguing that 
the main point of interest is the temporal behaviour of auto-correlation func- 
tions (for a detailed analysis of relation between auto-correlation functions 
and the fidelity criterion see [11]). 

To information encoded in an initial state p we associate a complex ob- 
servable A through 

p = p^(l + A) with Tr{ppA) = 0. (15) 
The observable A^ can then be used to recover this information as 

Tr(pA^) = Tr{ppAA^) > 0. (16) 
This detection, however, disappears in the long run 

hm Tr(^pe*^(At)) = Tr(p/3^^) = 0. (17) 
By the detailed balance property of C we have 

'Iv(pe*^(At)) =Tr(p^Ae*^(At)^ ^Tr(^^^e*^/2(A)e*(^*+^/2)(At)). (18) 
We then apply Schwarz's inequality to obtain 

Tr(p^Ae*^(At))|' = |Tr(p^e*^/2^A)e*(^'^+^/2)(At))|' (19) 
< {Tr(p^e*^/2(A)e*^/2(At)^|' (20) 

= {Tr(p^Ae*^(At))}'. (21) 
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Therefore the time auto- correlation functions of observables evolving under 
the dissipative part of the dynamics alone determine the useful lifetime of 
observables. As a detailed balance generator C is normal, the lifetime is 
determined by the smallest eigenvalue of —C different from 0. 

Generally, let 7i be a subspace of observables containing 1 and globally 
invariant under an ergodic detailed balance generator £, then 

Gap(^l^) 

:— min^|— A : 7^ A eigenvalue of j (22) 

= min(^|-(X, : X&n, ||X||/3 = 1, and (l,X)^ = o}^. 

(23) 



4 Master Hamiltonians 



For explicit computations it is convenient to rewrite the Davies generator 
in the standard Hilbert-Schmidt space. E.g., when dealing with composite 
systems, the Liouville scalar product introduces correlations between the 
different parties while the simplicity of the tensor structure is clearly visible 
in the Hilbert-Schmidt picture. Associating to an observable X the vector 

^■.^Xp% (24) 
we pass from Liouville to Hilbert-Schmidt space 

((^, ^) := Tr(^V- (25) 
The map (24) is unitary as (X , X)^ = ((/? , (/?). 

We can now unitarily transport the action of —Cauj, see (13), on Liouville 
space to an action kau, on Hilbert-Schmidt space 

Ku, = {Saiuj)L -VSaiuj)R}*{Saiuj)L - 5'c,(cj)^} 

+ {Sa{uj)R - T] Sa{uj)L}* {Sa{uj)R -r]Sa{uj)L} (26) 
= {{S„{uy)^ - T] {Saiujy)j^}{Saiuj)L -r]Sai^^)R} 

+ { (5«(a;)% - 7] {Sa{u;y)^}{Sa{u;)R -vSa{u^)L}. (27) 
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Here we introduced left and right multiplication 



if := X if and Xji ip := tpX 



(28) 



and 



T] := exp(— /3a;/2). 



(29) 



A * was also used instead of a f to distinguish between the Hermitian con- 
jugates for the usual and for the Hilbert-Schmidt inner product. 

By (26) each kau> is manifestly positive and the same is true for the sum 



The operator K is minus the unitary transform of the Davies generator jC 
under the map (24). In particular K and —C have the same eigenvalues, 
taking multiplicities into account. 

We have now reached a Hamiltonian picture of C in terms of K which 
is a sum of positive contributions with the remarkable property that these 
contributions have a common zero energy vector, the identity operator seen 

1 /2 

as the vector p^J in Hilbert-Schmidt space. This is rather unusual as the 
different terms in (30) generally don't commute, implying that X is a truly 
quantum Hamiltonian. This fictitious Hamiltonian, which will always be 
denoted by K, should be distinguished from the system Hamiltonian H^^^, 
we will call it master Hamiltonian. The positive constants g are not very 
relevant, we can in fact replace the g by arbitrary positive constants without 
changing the ground state of K. Even more, the full spectrum of K can be 
sandwiched between spectra of K^s with modified constants. 

As argued at the end of Section 3, the important parameter characterizing 
the relaxation time of the system is the spectral gap of the Davies generator 
C, see (22), which is equal to the ground state gap of the master Hamiltonian 
K as both are unitarily equivalent. Consider a non- negative matrix A with a 
non-trivial kernel. The smallest strictly positive eigenvalue g oi A is its gap 
above the kernel and will be shortly denoted by Gap( A) . We shall use in the 
sequel a number of simple facts about gaps. 

Lemma 1. i) Let A be a positive operator with non-trivial kernel, then any 

real number g such that A? > gA is a lower bound for Gap(A). 

a) Let A and B be positive operators such that Kqi[A-\- B) is non-trivial and 




(30) 



a u}>0 
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that Ker(A + B) = Ker(S), then Gap(A + B) > Gap(5). 

in) Let A and B be commuting positive operators such that KeT[A + B) is 
non-trivial, then Gap{A + B) > min(^{Gap{A), Gap(-B)}). 

We shall also need the following lemma: 

Lemma 2. For positive operators A and B, let A have gap g^ and {(p , B(j)) > 
gs for all normalized ip G Ker^A), then 



A + B> 



QaQb 



9A+\\B\ 

Proof. For any normalized vector ^0 we write 

ip = aip + hip' 



(31) 



(32) 



where 99 G Ker(A), G Ker(>l)-'- are normalized and \a\^ + \h\^ — 1. We 
then have 

(z^, {A + B)ij) > {r], Mr])>X_, (33) 
where (• , •) is the scalar product in C^, r] = (a, b), 



M 



"0 


0" 









+ 


9A_ 





{cp^Bcp) {ip.Bip') 
{ip',Bip) {ip',Bip') 



(34) 



and A_ is the smaller eigenvalue of the positive matrix M. The required 
estimate A_ is obtained from Lemma 3 below, noting that \\B\\ > \\M\\. □ 



Lemma 3. Consider the 2x2 matrix C given by 



y X 

X* z + u 



where u> Q and C is positive. Here, 

C 



(35) 



y X 
X z 



Let e_ be the smaller eigenvalue of C , then 

yu 



e_ > 



u+\\C'\ 



(36) 



(37) 
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5 The Ising ferromagnet on a ring 



The aim of this section is to analyze the relaxation to equilibrium of a system 
of Ising spins on a ring weakly coupled to a thermal environment. The 
model is too simple to be a reasonable candidate for a memory but it will 
provide useful techniques to deal with the 2D Kitaev model. 



5.1 The model 



The elementary building blocks of the model are two-level systems with the 
usual Pauli matrices and the identity as a basis of the observables 



(7^ = 



"0 1" 




"0 -I 


1 




i 



and cr^ = 



1 
-1 



(38) 



There are such spins arranged in a ring configuration, a subscript j = 
1,2, ... ,N denotes the position of the spin and + 1 is identified with 1. 

The interaction for an Ising ferromagnetic ring is given by bond observ- 
ables 

Zb := (7j(7|+i, with b = {j, j + 1}. (39) 
Clearly these bonds are not independent as they satisfy the relation 



1. 



(40) 



The ferromagnetic Ising Hamiltonian is then 

Rising _ J>0 



(41) 



where the prime on the summation reminds that we must take cyclic bound- 
ary conditions (40) into account. This spin model has a two-fold degen- 
erate ground state. The Zb, and any product of Z^s, all have expectation 1 
which determines the state completely but for one freedom in C^. Precisely 
this freedom is used to encode a single logical qubit. As before, the term 
qubit will be used to distinguish between the physical spins and an encoded 
abstract qubit which is exhibited by the following construction. 
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The observables 



X := (T^a^ • • • (T^ and Z := a'^ (42) 

commute with aU the bond observables and are therefore constants of the 
motion for the Hamiltonian evolution. Moreover, they satisfy the qubit rela- 
tions 

X = Z = Z\ X^ = = 1 and XZ + ZX = 0. (43) 

Actually, the Hamiltonian (41) could be taken more generically without al- 
tering the picture. Instead of (41) we could consider 

Rising J6>Cst. (44) 

b 

We would still have the same ground state and the constants of the motion 
would now be the algebra generated by the bond observables. The qubits (42) 
would still be preserved under this more generic kind of evolution. 

Let Q the be the qubit algebra generated by X and Z. We shall decompose 
the observables with respect to Q. An admissible bond on a ring of N 
sites is a configuration |b) = 62, • • • , ^at) with bj G {+, — } containing 
an even number of minuses. Next consider a Hilbert space Ti+ spanned 
by an orthonormal family {|b)} with b admissible. The algebra of linear 
transformations of T-C+ will be called the full bond algebra and denoted by 

"^rb ■ 

Lemma 4. The algebra of observables of the Ising model on a ring can be 
decomposed into the tensor product 

Q®X""- (45) 

Proof. A natural orthonormal basis of the Hilbert space of the full system 
consists of the tensor basis |e) where |e) = |ei, £2, . . . , ctv) corresponds to the 
eigenvalues of all cr^'s. An alternative labelling is to specify the eigenvalue of 
the qubit Z and |b) of all the admissible bond observables. A^^^ is then gen- 
erated by the bond observables and by all operators that flip two intersecting 
bonds and leave Q untouched, i.e. by crj with j 7^ 2. □ 
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5.2 The gap of the generator 



Let us now turn to the ring in a thermal environment. A generic interaction 
Hamiltonian is of the form 

N N N 

//-t = J] ® /, + J] ® /^+, + J] a| ® hN+,. (46) 
j=i j=i j=i 

where fj is a self-adjoint field of the j-th bath. To each of the terms in (46) 
there corresponds a number of terms in the Davies generator labelled by the 
Bohr frequencies of the coupling operators. As explained in Section 2 every 
contribution is a positive operator on Liouville space and the identity is a 
common element of their kernels. 

As H^^^, see (41), is a sum of local terms which commute amongst them- 
selves, the Hamiltonian system dynamics inherits this strong locality prop- 
erty: an observable living at site j will at most spread to the sites j — 1, j, 
and j + 1- We explicitly compute the evolution of the crj coupling operators 

The Fourier components of crj are given in terms of projection operators 

P° = ^ (1 - Z,Zk>) and P± = I (1 T Z,) (1 T Zk>) (48) 
with j = {6, b'}: 

Note that the projectors select neighbouring bond configurations that are 
equal (up or down) or different and that they satisfy P° -|- P+ -|- P" — 1. 
The Bohr frequencies of the Ising Hamiltonian are 0, ±4 J, ±8 J, . . . but only 
the Bohr frequencies and ±4J contribute to the evolution of the coupling 
operators crj, this is due to the strong locality property of the Ising dynam- 
ics. Similar expressions hold for the crj terms, it suffices to replace x by 
y in (47-49). As the detailed structure of the bath operators fj is rather 
irrelevant and to avoid inessential constants we shall choose the ^^(a;) in 
the Davies generator (6) in a convenient way, always assuming that they are 
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strictly larger than a positive constant independent of the system size and 
the frequency. To be specific, the constants are chosen in such a way that 
the contributions coming from crj (or cr|) are of the following form: 



£^(A) = i a] [A , ttj] + h+ [a] , A] + h_ aj [A , a]] 
+ h^[aj,A]a]}-^ho[al[alA]]. 



(50) 



so that hj^ + h_ = 2ho. Here 

7 = (52) 

Theorem 1 (Gap Ising ring). Let L he the Davies generator for a system 
of N Ising spins on a ring, see (41), obtained by the weak-coupling procedure 
from the interaction Hamiltonian i?'"* as in (46), then 

Gap(£) > \ e-^^-^. (53) 

Proof. In Section 4 we showed that the gap of the Davies generator coincides 
with the gap of the master Hamiltonian which consists of a sum of positive 
contributions with a unique common zero-energy eigenvector 1 in Liouville 
space. By Lemma 1 ii) we will only lower the gap by dropping terms in the 
master Hamiltonian (30) provided we don't introduce additional zero-energy 
vectors in doing so. We shall use this property to obtain a simpler expression 
for a lower bound on the gap. In fact we shall only retain terms coming 
from (Tj , j = 1,2, . . . , N and af. This leads to the following structure for the 
simplified generator 

C^,+C\ + L (54) 

Here refers to the contribution coming from a^, C\ to that from a\ and 
C from the contributions of the with j = 2,3, . . . , N. 

Clearly, by (42) all crj with j = 2, 3, . . . , iV belong to and so do the 
projectors P° and for all j. Therefore the Davies operators aj and a° with 
j 1 belong to the algebra .4^"^^ This imphes that C leaves Q untouched 
and that A^^^ is globally invariant. By abuse of notation we shall from now 
on use the notation id instead of C 
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Next we consider the orthogonal decomposition of Q into a direct sum 
of four ID spaces, generated by 1, X, Y and Z respectively. Here Y = iZX. 
Note that this is an orthogonal decomposition not only for the Liouville 
scalar product but for the Hilbert-Schmidt scalar product as well. This 
decomposition induces one of the full algebra of observables as 

(i ® X"'') © (x ® X"") ® (y ® X'') © (z ® X"'') • (55) 

Clearly id^C leaves each of these blocks invariant and has a same action 
on the X"^^ factors. Moreover, both Cf and C\ are block-diagonal in the 
decomposition (55). This can be seen by considering their explicit actions 
on elements of the type Ai^ B with A e {1, X, Y, Z} and B e X"^^- 

As a first simplification, we observe that each of the three terms in (54) 
is negative on X A^^\ Y (g) X"^^ and Z (g) X"""- We further simphfy our 
description by retaining only on X (8) X"^^ >Cf on Y (8) X"" either one 
of the two on Z (g) X"^^- 

We shall now bound from below —C — Cf restricted to Z (g) X"^^ the two 
other parts can be handled in the same way. By Proposition 1, which will be 
proved in Section 5.3, C is ergodic and gapped with the estimate given 
by (63). Hence C restricted to Z ® X"^^ has gap bounded by /i_/2 above its 
kernel Z (g) 1. Lemma 2 then says that on this subspace 

-C-n> ^ \ — -. 56 

Noting that 

G := h+P^ + h-P^ + hoPi > /i- (57) 
with h±, ho given by (51), we obtain 

-(Z® 1, £^ Z® 1) =2(1, G) > 2/i_ and \\Cl\\<2. (58) 
This gives us the following estimate 

(-1 - £1) > y-^. (59) 

For the first term 1 ® X"" (55) we retain only id <S>C As said, £ is 
ergodic and gapped with the bound /i_/2. We therefore conclude that 

Gap(£) > -n({-, ^}) - ^ > - > ^ (60) 

2 ' 2 ' 
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since 1 > /i- > 7^. □ 



5.3 A bond model 

In this section we consider the spectral gap of the generator L introduced in 
the proof of Theorem 1. As before, let be the Hilbert space spanned by 
the basis vectors |b) labelled by bonds on a ring, i.e. |b) = |6i,62, • • • , &iv) 
with hj e {+, — } and |b) contains an even number of minuses. Consider the 
following operators on 7^+ 

«i = (!--)(+ + 1)66', «l = (l + +)(--|)66', and 

(61) 



«;=(i-+)(+-iw + (i + -)(-+i) 



w 



where the bonds h and h' satisfy j — {6, 6'}. We are interested in the spectral 
gap of the operator L acting on the linear transformations B{'H'^) of 

1 ^ 

^ j=2 (62) 
Proposition 1. T/ie generator C is ergodic and gapped with bound 



Gap(^) > T = T^V ^^^^ 



and 7 defined in (52). 



In order to estimate the gap of C we shall use details on its structure as an 
operator on Liouville or Hilbert-Schmidt space. For any A C {1,2,..., A^} 
with (A^ — #(A)) even we consider the subspace i3(7i+)(A) of i3(7i_)_) spanned 
by the rank one operators |b)(b'| where b and b' coincide on A and where 
all the bonds of b' sitting on the complement of A are flipped with respect 
to these of b. These subspaces are orthogonal, both with respect to the 
Liouville and the Hilbert-Schmidt scalar products, so that we can write 

B{n+)^®B{n+){k). (64) 

A 
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It is furthermore not hard to see that 

B{n+){K) = ^"^(A) ® J^(A). (65) 

Here ^'^'^(A) is the AbeUan algebra generated by the projectors on bonds 
restricted to A and J^(A) is the space spanned by proper flips of bonds re- 
stricted to the complement of A. It is straightforward to check that every 
block B{T-C+){A) is left invariant by £,. Due to the positivity of JC we obtain 
the following structure for the generator. 

Lemma 5. The operator C, see (62) restricted to the algebra B{H+) is block 
diagonal with respect to the decomposition (64)- 

Proof of Proposition 1. We rely on the block diagonal structure of C as de- 
scribed in Lemma 5. The proof consists of the following steps 

1) obtaining a lower bound for minus the restriction of C to blocks with non- 
trivial flip part, which follows immediately this enumeration 

2) obtaining a lower bound for the gap of the generator >Cab which is the 
restriction of C to the subspace with trivial flip part, see Proposition 2 

3) proving the ergodicity of £, see Lemma 6. 

A partition of {1, 2, ... , N} into A and its complement induces a partition 
of the spins j — {6, 6'} into three sets: Fflip, Fab and Tint, consisting of pairs 
of neighbouring bonds from A*^, pairs of bonds from A and pairs that have 
one bond in A and the other in its complement. This partition induces 
accordingly a decomposition of the generator C. By explicit use of (61) one 
verifles that for j e Fab the flip part is left untouched while the Abehan 
part is mapped into itself. The same holds true for j e Fflip. It remains 
to consider contributions coming from j G Fint- Wc just give an explicit 
example. Suppose that j = {b.b'} with 6 G A and b' G A'^ such as (| -|- 
+)(+- Dw- Using (61) weflnd 

++)(+-!) = —(h. + ho)\ ++){+-]. (66) 

with h-, ho given by (51). It is then straightforward to verify that 

->C > #(Fint)/i- > /i-. (67) 

We should therefore distinguish three cases: the Abehan case Fflip = 0, see 
Proposition 2, the pure flip case Fab = 0, see the remarks following the proof 
of Proposition 2, and the case Pint 7^ which was dealt with above. □ 
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Let be the restriction of the generator C to the diagonal algebra 
generated by projectors of the form |b)(b| with b an admissible bond on the 
ring. As the Gibbs density matrix pjs is also diagonal these projectors are 
an orthogonal family. Minus the restricted generator C^h can now easily be 
expressed with respect to the Hilbert-Schmidt inner product 



K.. 



ab 



(68) 



hV 



where 



hb' = 











V- 



1+7^ 







1 

2 





1+^2 \ 




1 



(69) 



1+7^ 



/ 



Due to our choice of constants Qaif^j) (as done in 50) are projectors. The 
summation extends over an open chain of bonds, indeed the contribution of 
the site j — 1 oi the ring was removed. 



Proposition 2 (Gap on Abelian factor). With the notation of above 



Gap(£ab) = Gap(Xab) > - 



r 



+ 7' 



(70) 



Proof. A lower bound for the gap of K^h can be obtained by comparing 



(Xab) with a multiple of K^h-, see Lemma 1 i). We start out by writing 

(-^ab)^ = (^12 + /i;23 H V kN-lN^ 



— + ^23 + • • ■ + k'^N-lN + ^12^23 + ^23^12 + 

= A;i2 + A;23 H \- kw-iN + kuhs + k2zki2 + 



(71) 

(72) 
(73) 



where we used that kw is a projector. Suppose now that we find a positive 
g such that 

ki2 + A;23 + 2A;i2/c23 + 2/c23A;i2 > g{ki2 + A;23) , (74) 

then we may obtain a lower bound for K^^. Projectors k that have no bond 
in common commute, implying that products of such projectors are positive. 
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We may also decrease the coefficient of a projector. This leads to 

(i^ab)' > ku + k23 + --- + ki2k23 + k23ku + ■■■ (75) 

> ki2 + ^k23 + ki2k23 + k22.kl2^ 

+ (^^ ^23 + ^ ^34 + k23ku + k2Ak23^ H 

+ ^2 ^N-2N-\ + 2 kN-lN + kN-2N-lkN-lN + ArA;Ar-2 AT^ 

(76) 

> { (2 ^12 + 2 ^23) + (2 ^23 + - ^34) H 

) } (77) 

>f^ab. (78) 

It remains to obtain a g satisfying (74). As the /c's are projectors, the con- 
dition is equivalent to 

{k,2 + k23Y>^{ku + k23). (79) 

The eigenvalues of ki2 + /C23 are {0, 2, 2(1+^2^ > 2(1+7^) }• '^^^^ leads to the 
estimate 

Gap(4b) > (80) 

□ 

It remains to investigate the case Fab = 0, i-e. to consider the action of 
jC on the space of linear combinations of matrix units of the form |b)(b'| 
where b is an admissible configuration of bonds and where b' is the bond 
configuration obtained from b by flipping every bond. The configuration 
b' should also be admissible which is impossible for odd N. For even N 
essentially a similar approach to the one in Proposition 2 can be used to 
obtain that —C > 1/2 on this subspace. In order to keep the proof concise 
we shall assume in the sequel that is odd. Finally, to finish the proof of 
Proposition 1, it remains to show the ergodicity of C. 



Lemma 6 (Ergodicity). The generator C is ergodic on B{7i 



+) 
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Proof. As in the analysis of the gap of the generator, we shall assume that 
is odd, this simplifies the arguments. The general case can however be 
handled as well. The commutant of the Davies operators a° and a^, j = 
2, 3, . . . , A" is the same as the commutant of {crj, Zj, + Zf,' with j = {b, 6'}}. 
As A^ is odd this turns out to be the commutant of B{H+) which is trivial. □ 

6 Kitaev's 2D model 

Even if Kitaev's 2D model is more complicated than the Ising model on a 
chain many ideas in the analysis of it's thermal stability are quite similar, cer- 
tainly on the mathematical level. We shall stress this similarity by following 
the approach used for the Ising case. 

6.1 The model 

Kitaev's 2D model lives on a L x L toroidal lattice A. The microscopic 
constituents are spins located at the midpoint of the edges of the lattice and 
represented by filled and open dots in Fig. 1, so there are 2L^ of them. The 
interactions between the spins are given by star and plaquette terms Xg and 
Zp. A star is a cross whose vertical vertices lie on the lattice of black dots, 
while the vertical vertices of a plaquette lie on the lattice of white dots. The 
actual star and plaquette observables are 

= n Zp^Yl a]. (81) 

Stars and plaquettes correspond to the grey shapes in the figure. As stars and 
plaquettes have either or 2 sites in common, [X^, Zp] = 0. So, the algebras 
«4.x, ^z, and ^xz generated by the X^, the Zp, and by both together are 
Abelian. 

The Hamiltonian of the model is 

i/ft = JX, - ^ JZp, J>0. (82) 

Similarly to the Ising model the ground states are totally unfrustrated: all 
Xg and Zp have expectation 1. This is actually not sufficient to fully de- 
termine the state of all spins as the star and plaquette observables are not 
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Figure 1: Kitaev's lattice 



independent: because of the periodic boundary conditions they satisfy 

= 1 and Y[Zp = 1. (83) 

s p 

As a consequence two topological qubit freedoms are left which may be used 
for encoding. The Hamiltonian (82) can be chosen more generically by multi- 
plying the individual star and plaquette observables by positive but otherwise 
arbitrary coefficients, this will not change the set of ground states. Here too, 
it is natural to consider the commutant of such a generic Hamiltonian which 
consists of a product of two qubit algebras and Axz- This is seen quite ex- 
plicitly by introducing, similarly to (42), observables for two encoded qubits 



(84) 



Here Ci, rfi, C2 and c?2 are the loops shown in Fig. 1. Unlike for the Ising 
ring, all qubit observables are very delocalized. 

Let us divide the set of spins into four disjoint subsets, sec Fig. 2: the 
snake, the comb, spin 1 and spin 2. Spin 1 is located at the crossing of Xi 
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Figure 2: Partitioning Kitaev's lattice 



and Zi and similarly for spin 2. Note that the qubit Xi has been modified a 
little, so that it closely follows the snake. 

Let ^p^" be the algebra is generated by all plaquettc obscrvables and 
by the crj where j belongs to the snake while A^^^^ is generated by the star 
observables and by the (7| with j belonging to the comb. It is obvious from the 
construction that these algebras commute and that they are both isomorphic 
to a full matrix algebra of dimension 2^ Indeed, considering A^^^^ first 
we have — 1 independent commuting plaquette observables. Admissible 
fiips of these plaquette observables should comply with condition (81). Now 
any such flip can be obtained by concatenating flips induced by a| with j 
belonging to the snake. A similar reasoning holds for A^^\ Therefore the 
commutant of these two algebras is a 4 dimensional matrix which may be 
identifled with the tensor product of two qubit algebras Qi and Q2 generated 
by the logical qubits (84). We have therefore shown that 

Lemma 7. The algebra of observables of the 2D Kitaev model on a toric 
L X L lattice can be decomposed into the tensor product 

Qi (g) Q2 ® ^p"" ® A'f. (85) 

There is a striking analogy between Lemma 7 and Lemma 4. The plaque- 
ttes or stars should be compared to the bonds. Moreover, the operators that 
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flip plaquettes or stars are located on essentially one-dimensional structures: 
the snake and the comb. This makes the 2D Kitaev model essentially iso- 
morphic to two Ising like models on ID structures. For the factors, A^^^^ and 
^g"" we shall therefore consider the same block decomposition as for A^^^: 

A'f = © AfiA) J^p(A) and = ^'''(A) ® ■^s(A). (86) 

Moreover, the factor ^p"" is fully isomorphic to the full bond algebra A^^^ of 
the ID Ising model, hence we shall further on use the same notation, keeping 
in mind that bonds translate into plaquettes. 



7 Gap for Kitaev's model 

In the Davies generator we now retain only contributions coming from crj 
and (t|. The generator has the form 

N N 

C = Cxj + Cz,j = Cx + Cz- (87) 

where j runs over all qubits. Note that the action of C is ergodic and that 
Cx and Cz commute. Therefore it is enough to show that e.g. Cx has a gap. 

We first note that Cx is block diagonal, where the blocks are even finer 
than those determined by formulas (86). Namely, we consider for some fixed 
A blocks of the form As0 A^p^^ (8 Z,^X^, where is a proper flip from 
J^s- One then finds, similarly to formulas (91)-(94) in [4]), that 

Cx{As ® ® A'f ® Z,X^) = A ® i^s ® ^'xiA'^'') ® Z,X^ 



Here C'^ is defined as follows: for any j which belongs to the symmetric 
difference of S, i.e. the set of qubits where Fg has ai matrices, and the 
support of Zu the term C'^ j{A) has the form — ^ (a^aA -|- Aa'^a) — a'^Aa. 
This is so because for such a j, the cr^ term from the generator has to be 
commuted over a cr^ occurring in Fg or in Zj,. This results in flipping the sign 
of the contribution a^{A)a. For all other j's, the generator is unchanged: 
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We now have to show that has a gap on each block. To this end we 
introduce which consists solely of contributions from sites belonging 

to the snake, see Fig. 2. We shall use that 

L^:^''' (A ®F,® ® Z,X^) = A i^s C^T'^' (4"") ® Z,X^ (89) 
and that 

Gap(£r'^)|^^^^^>^ and Ker(£r'^) |^^^^_ = {1}. (90) 

Here /i_ is the same constant as in (51). This follows from Theorem 1 since 
^snake restricted to ^p"^^ is fully isomorphic the Ising model on 

We shall now consider three types of blocks. For the first two we shall 
bound —Lx from below, exactly like in the Ising model for blocks of type 
Z®^^"". In the third type of block —Lx will be gapped by the same estimate 
as in (70). 

Case i). Consider all blocks where Zi is present, the case where Z2 is 
present is analogous. We then remove lot of interactions and retain 

4' = £r'^+£^^. (91) 

Note that has support outside the comb. We have therefore 
C'^{A, (g) (g) ZiX^) =A ® i^s ® >Cr''(4"'^) ® ZiX^ 

We now apply Lemma 2. From discussion below (88) it follows that 

(Zi,£,,(Zi)) > 2(1,(7) ;3>2/i_ (93) 

with G as in (57). Within the block under consideration we obtain exactly 
the same bound as for the Ising model in Theorem 1 

-4' > y (94) 

where we have used ||LgJ| < 2. This implies the same estimate for the full 
generator on this block 

-Cx > y. (95) 
23 



(92) 



Case a). Consider now blocks where Z^, is not present but Fs is. We then 
take 

Cl = CT^^ + (96) 

where go is any spin belonging to the set S determined by the flip Fg. A 
similar argument as for Case i) then yields the same estimate (95). 

Case in). Finally, consider blocks where neither F^ nor Z^, is present. 
We then simply take C^^^^^^ itself. For this kind of block C^^'^^^ has a gap ^ 
above its kernel which is generated by As ® X^. One immediately sees that 
this coincides with the kernel of jCx on this block. Therefore by Lemma 1 (ii) 
we have 

Gap(>C^) > ^. (97) 

Combining the three cases, we conclude that Cx has a gap above its kernel 
which satisfies ^ 

Gap(£,) > ^. (98) 

We have therefore completed the proof of the following theorem: 

Theorem 2 (Gap Kitaev's 2D model). Let C be the Davies generator for 
the 2D Kitaev model of size L x L, then 

Gap(£) > \ e"^^^. (99) 



8 Conclusion 



We presented in this paper a rigorous analysis confirming the heuristic ar- 
guments in [2, 3] that the 2D Kitaev model is not stable with respect to 
thermal fluctuations. The physical mechanism behind this instability is the 
absence of interactions between the thermal excitations. Indeed, Kitaev's 2D 
model can be translated into a non interacting quasi-particle model (anyons), 
similarly to non-interacting kinks in the ID Ising model. Therefore, except 
for a gas of bosons where the statistics plays a major role, only systems with 
truly interacting thermal particles like the 2D Ising or the 3D and 4D Ki- 
taev models could support delocalized metastable observables. In this paper 
the analysis of the dissipative generator is translated into the analysis of the 
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ground state of a quantum spin Hamiltonian K called master Hamiltonian. 
Quantum spin chain techniques can then be used to obtain lower bounds on 
the spectral gap of K, leading to instability results as in this paper. 

The same technique can however also be applied to obtain for suitable 
models upper bounds for the spectral gap of K and to unravel the form of 
possible metastable observables. Such models and their relevance for the im- 
plementation of quantum memory is the topic of a forthcoming publication. 
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